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We give a comprehensive presentation of methods for calculating the Casimir 
force to arbitrary accuracy, for any number of objects, arbitrary shapes, suscepti- 
bility functions, and separations. The technique is applicable to objects immersed 
in media other than vacuum, nonzero temperatures, and spatial arrangements in 
which one object is enclosed in another. Our method combines each object's classi- 
cal electromagnetic scattering amplitude with universal translation matrices, which 
convert between the bases used to calculate scattering for each object, but are other- 
wise independent of the details of the individual objects. The method is illustrated 
by re-deriving the Lifshitz formula for infinite half spaces, by demonstrating the 
Casimir-Polder to van der Waals cross-over, and by computing the Casimir interac- 
tion energy of two infinite, parallel, perfect metal cylinders either inside or outside 
one another. Furthermore, it is used to obtain new results, namely the Casimir 
energies of a sphere or a cylinder opposite a plate, all with finite permittivity and 
permeability, to leading order at large separation. 



I. INTRODUCTION 

Materials that couple to the electromagnetic field alter the spectrum of its quantum and 
thermal fluctuations. The resulting change in energy depends on the relative positions of 
the objects, leading to a fluctuation-induced force, usually called the Casimir force. This 
force has been the subject of precision experimental measurements [UlSlElllllSliniElIHlini 
[TU| [TT| \T2\ [T5| dH [T51 UHl [17] and can influence the operation of nanoscale devices dS] . 

Casimir and Polder calculated the fluctuation-induced force on a polarizable atom in front 
of a perfectly conducting plate and between two polarizable atoms, both to leading order at 
large separation, and obtained a simple result depending only on the atoms' static polariz- 
abilities [19]. Casimir then extended this result to his famous calculation of the pressure on 
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two perfectly conducting parallel plates [20]. Feinberg and Sucher |2T1 122] generalized the 
result of Casimir and Polder to include both electric and magnetic polarizabilities. Lifshitz, 
Dzyaloshinskii, and Pitaevskii extended Casimir's result for parallel plates by incorporating 
nonzero temperature, permittivity, and permeability into a general formula for the pressure 
on two infinite half-spaces separated by a gap [231 IMl ESI U7\ ■ 

In order to study Casimir forces in more general geometries, it turns out to be advanta- 
geous to describe the influence of an arrangement of objects on the electromagnetic field by 
the way they scatter electromagnetic waves. In this article we derive and apply a represen- 
tation of the Casimir energy, developed in Refs. [281 129], that characterizes each object by 
its on-shell electromagnetic scattering amplitude. The separations and orientations of the 
objects are encoded in universal translation matrices, which describe how a solution to the 
source-free Maxwell's equations in the basis appropriate to one object looks when expanded 
in the basis appropriate to another. The translation matrices depend on the displacement 
and orientation of coordinate systems, but not on the nature of the objects themselves. The 
scattering amplitudes and translation matrices are then combined in a simple algorithm 
that allows efficient numerical and, in some cases, analytical calculations of Casimir forces 
and torques for a wide variety of geometries, materials, and external conditions. We will 
generalize the formalism summarized in [28] further to show how it applies in a wide variety 
of circumstances, including: 

• n arbitrarily shaped objects, whose surfaces may be smooth or rough or may include 
edges and cusps; 

• objects with arbitrary linear electromagnetic response, including frequency-dependent, 
lossy electric permittivity and magnetic permeability tensors; 

• objects separated by vacuum or by a medium with uniform, frequency- dependent 
isotropic permittivity and permeability; 

• zero or nonzero temperature; 

• and objects outside of one another or enclosed in each other. 

These ideas build on a range of previous related work, an inevitably incomplete subset of 
which is briefly reviewed here: Scattering theory methods were first applied to the parallel 
plate geometry, when Kats reformulated Lifshitz theory in terms of refiection coefficients 
|30j . Jaekel and Reynaud derived the Lifshitz formula using reflection coefficients for lossless 
infinite plates [31] and Genet, Lambrecht, and Reynaud extended this analysis to the lossy 
case [3^. Lambrecht, Maia Neto, and Reynaud generalized these results to include non- 
specular reflection [33] . 

Around the same time as Kats's work, Balian and Duplantier developed a multiple scat- 
tering approach to the Casimir energy for perfect metal objects and used it to compute 
the Casimir energy at asymptotically large separations [311 ES] at both zero and nonzero 
temperature. In their approach, information about the conductors is encoded in a local 
surface scattering kernel, whose relation to more conventional scattering formalisms is not 
transparent, and their approach was not pursued further at the time. One can flnd multiple 
scattering formulas in an even earlier article by Renne [36], but scattering is not explicitly 
mentioned, and the technique is only used to rederive older results. 

Another scattering-based approach has been to express the Casimir energy as an integral 
over the density of states of the fluctuating fleld, using the Krein formula [371 EH El] to 
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relate the density of states to the iS-matrix for scattering from the ensemble of objects. This 
iS-matrix is difficult to compute in general. In studying many-body scattering, Henseler and 
Wirzba connected the iS-matrix of a collection of spheres POJ or disks [41j to the objects' 
individual iS-matrices, which are easy to find. Bulgac, Magierski, and Wirzba combined 
this result with the Krein formula to investigate the scalar and fermionic Casimir effect for 
disks and spheres ^21 H3l HI] . Casimir energies of solitons in renormalizable quantum field 
theories have been computed using scattering theory techniques that combine analytic and 
numerical methods j45j. 

Bordag, Robaschik, Scharnhorst, and Wieczorek P6| Wl\ introduced path integral meth- 
ods to the study of Casimir effects and used them to investigate the electromagnetic Casimir 
effect for two parallel perfect metal plates. Li and Kardar used similar methods to study 
the scalar thermal Casimir effect for Dirichlet, Neumann, and mixed boundary conditions 
llHllin]- The quantum extension was developed further by Golestanian and Kardar [501 \5T\ 
and was subsequently applied to the quantum electromagnetic Casimir effect by Emig, 
Hanke, Golestanian, and Kardar, who studied the Casimir interaction between plates with 
roughness |^ and between deformed plates [53] . (Techniques developed to study the scalar 
Casimir effect can be applied to the electromagnetic case for perfect metals with translation 
symmetry in one spatial direction, since then the electromagnetic problem decomposes into 
two scalar ones.) Finally, the path integral approach was connected to scattering theory by 
Emig and Buescher [54] . 

Closely related to the work we present here is that of Kenneth and Klich, who expressed 
the data required to characterize Casimir fluctuations in terms of the transition T-operator 
for scattering of the fluctuating field from the objects [55] . Their abstract representation 
made it possible to prove general properties of the sign of the Casimir force. In Refs. 
[28] [29] , we developed a framework in which this abstract result can be applied to concrete 
calculations. In this approach, the T-operator is related to the scattering amplitude for each 
object individually, which in turn is expressed in an appropriate basis of multipoles. While 
the T-operator is in general "off-shell," meaning it has matrix elements between states with 
different frequencies, the scattering amplitudes are the "on-shell" matrix elements of this 
operator between states of equal frequency.^ In this approach, the objects can have any shape 
or material properties, as long as the scattering amplitude can be computed in a multipole 
expansion (or measured). The approach can be regarded as a concrete implementation of 
the proposal emphasized by Schwinger [56] that the fluctuations of the electromagnetic fleld 
can be traced back to charge and current fluctuations on the objects. This formalism has 
been applied and extended in a number of Casimir calculations [STJ [581 EHl EQ] [611 [62] . 

The basis in which the scattering amplitude for each object is supplied is typically asso- 
ciated with a coordinate system appropriate to the object. Of course a plane, a cylinder, 
or a sphere would be described in Cartesian, cylindrical, or spherical coordinates, respec- 
tively. However, any compact object can be described, for example, in spherical coordinates, 
provided that the matrix of scattering amplitudes can be either calculated or measured in 
that coordinate system. There are a limited number of coordinate systems in which such a 
partial wave expansion is possible, namely those for which the vector Helmholtz equation 
is separable. The translation matrices for common separable coordinate systems, obtained 



Because of this relationship, these scattering amphtudes are also referred to as elements of the T-matrix. 
In standard conventions, however, the T-matrix differs from the matrix elements of the T-operator by a 
basis-dependent constant, so we will use the term "scattering amplitude" to avoid confusion. 
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from the free Green's function, are supplied in Appendix [Cj For typical cases, the final com- 
putation of the Casimir energy can be performed on a desktop computer for a wide range 
of separations. Asymptotic results at large separation can be obtained analytically. 

The primary limitation of the method is on the distance between objects, since the basis 
appropriate to a given object may become impractical as two objects approach. For small 
separations, sufficient accuracy can only be obtained if the calculation is taken to very high 
partial wave order. In the case of two spheres, the scattering amplitude is available in 
a spherical basis, but as the two spheres approach, the Casimir energy is dominated by 
waves near the point of closest approach [63]. As the spheres come into contact an infinite 
number of spherical waves are needed to capture the dominant contribution (see Section 



III for further discussion). A particular basis may also be fundamentally inappropriate at 



small separations. For instance, if the interaction of two elliptic cylinders is expressed in 
an ordinary cylindrical basis, when the elliptic cylinders are close enough one may not fit 
inside the smallest circular cylinder that encloses the other. In that case the cylindrical 
basis would not "resolve" the two objects (although an elliptic cylindrical basis would; see 
Section III). Finally, for a variety of conceptual and computational reasons, we are limited 
to linear electromagnetic response. 

To illustrate this general formulation, we provide some sample applications, including the 
closed-form expressions for computing the interaction of a plate and a sphere with finite, 
uniform, frequency-dependent electric permittivity and magnetic permeability. We present 
the Casimir interaction energy explicitly at asymptotically large separations in terms of the 
zero frequency permittivities and permeabilities of the objects. Although most experiments 
have centered around the sphere-plate configuration [H El [3l E], El El HOl [121 UHl EH EH], it 
is only recently that the force between a dielectric sphere and an idealized metallic plate 
has been obtained for all distances |66]. Subsequently, this result has been extended to the 
situation where both objects are described by the plasma model [67]. In addition, we present 
the Casimir interaction energy of a plate and a cylinder at asymptotically large distances in 
terms of the two objects' zero frequency permittivities and permeabilities. Results beyond 
the leading order using our closed-form formulation are not explicitly included, but all the 
essential formulas are contained here. These results extend the perfect metal cylinder and 
plate results presented in Ref. [HS]- 

The article is organized as follows: In Section II we review the derivation of the ground 
state energy of a field theory using path integrals. In Section III we expand the free electro- 
magnetic Green's functions in terms of regular and outgoing waves, taking into account that 
the pairs of waves in the expansion are evaluated with respect to two different coordinate 
systems. This analysis yields the translation matrices. Section IV provides an overview of 
elements of scattering theory we will use, including the connection between the T-operator 
and the scattering amplitude. In Section V the path integral expression for the energy is 
re-expressed in terms of the results of the preceding two sections, yielding the main result, 
Eq. (V.13). In Section VI sample applications are presented: A short derivation of the 
Lifshitz formula, the cross-over between van der Waals and Casimir regimes for two atoms, 
a general derivation of previous results for cylinders j69l |70] , and new results for the energy 
between a dielectric sphere or cylinder and a dielectric plane. A discussion follows in Section 
VII. 



II. CASIMIR ENERGY FROM FIELD THEORY 
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A. Electromagnetic Lagrangian and action 

We consider the Casimir effect for objects without free charges and currents but with 
nonzero electric and magnetic susceptibilities. The macroscopic electromagnetic Lagrangian 
density is 

/: = 1(E.D-B-H). (II.l) 

The electric field E(t,x) and the magnetic field B(t, x) are related to the fundamental 
four-vector potential A'^ by E = —c~^dtA — and B = V x A. We treat sta- 

tionary objects whose responses to the electric and magnetic fields are linear. For such 
materials, the D and B fields are related to the E and H fields by the convolutions 
D(t, x) = dt' e{t', x)E(t-t', x) and B{t, x) = dt' fi{t', x)H(t-t', x) in time. We con- 
sider local, isotropic permittivity and permeability, although our derivation can be adapted 
to apply to non-local and non-isotropic media simply by substituting the appropriate non- 
local and tensor permittivity and permeability functions. A more formal derivation of our 



starting point Eq. (II.l), which elucidates the causality properties of the permeability and 
permittivity response functions, is given in Appendix [A| 

We define the quantum-mechanical energy through the path integral, which sums all 
configurations of the electromagnetic fields constrained by periodic boundary conditions in 
time between and T. Outside of this time interval the fields are periodically continued. 
Substituting the Fourier expansions of the form E(t, x) = J2'^=-oo ^i^n, x)e~*'^"* with Un = 
27in/T, we obtain the action 

S{T) = - [ /" dx (E ■ D - B ■ H) = -T V /" rfx (E* ■ eE - B* ■ fi-^B) , (II.2) 

where e, E, /i, and B on the right-hand side are functions of position x and frequency Un, 
and we have used D(c<j,x) = e(co', x)E(co', x) and 11(0;, x) = x) -^(^'-^)- 

From the definition of the fields E and B in terms of the vector potential A^^, we have 
V X E = i-B, which enables us to eliminate B in the action. 



^ oo 
n=—oc 

where 



dx 



E* ■ ( I X Vx E ^E* ■ VE 



(11.3) 



V = I^(l-eK,x)) + Vx r-1 Vx (II.4) 

is the potential operator and we have restored the explicit frequency dependence of e and 
/i. The potential operator is nonzero only at those points in space where the objects are 
located (e 7^ 1 or 7^ 1). At small frequencies, typical materials have e > 1 and /i ~ 1, and 
V can be regarded as an attractive potential. 

In the functional integral we will sum over configurations of the field . This sum must 
be restricted by a choice of gauge, so that it does not include the infinitely redundant gauge 
orbits. We will choose to work in the gauge = 0, although of course no physical results 
depend on this choice. 
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B. Casimir energy of a quantum field 

We use standard tools to obtain a functional integral expression for the ground state 
energy of a quantum field in a fixed background described by V(a;, x). The overlap between 
the initial state |Ea) of a system with the state |Eb) after time T can be expressed as a 
functional integral with the fields fixed at the temporal boundaries [7T] , 

(E,|e-™|E„) = [VA |k,.o)^e. e^'^K (II.5) 



where S{T) is the action of Eq. (II. 2) with the time integrals taken between zero and T, and 
H is the corresponding Hamiltonian. 

If the initial and final states are set equal and summed over, the resulting functional 
integration defines the Minkowski space functional integral 

Z{T) = ^(Eje-^^^/^|E,) = tr e"^^'^/'^ = f VAe^'^K (II.6) 

which depends on the time T and the background potential Y{uj, x). The partition function 
that describes this system at temperature is defined by 

Z{f3) = Z{-ihf3) = tr e"^^, (II.7) 

and the free energy F of the field is 

F{P) = -^\ogZiP). (II.8) 

The limit /3 ^ oo projects the ground state energy out of the trace, 

So = F{/3 = oo) = - lim ]- log Z{(3). (II.9) 

The unrenormalized energy So generally depends on the ultraviolet cutoff, but cutoff- 
dependent contributions arise from the objects individually and do not depend on their 
separations or orientations. Such terms can also arise after renormalization if objects are 
assumed to constrain electromagnetic waves with arbitrarily high frequencies (for example, 
if the fields are forced to vanish on a surface). Such boundary conditions should be regarded 
as artificial idealizations; in reality, when the wavelengths of the electromagnetic waves be- 
come shorter than the length scales that characterize the interactions of the material, the 
influence of the material on the waves vanishes [72]. Accordingly, the potential V should 
vanish for real materials in the high-frequency limit. Since we are only interested in energy 
differences, we can remove these divergences by subtracting the ground state energy of the 
system when the objects are in some reference configuration. In most cases we will take 
this configuration to have the objects infinitely far apart, but when calculating Casimir en- 
ergies for one object inside another, some other configuration must be used. We denote the 
partition function for this reference configuration by Z. In this way we obtain the Casimir 
energy, 

£ = - lim ^log Z{[3)/Z{(3). (11.10) 

Throughout our calculation of we will thus be able to neglect any overall factors that are 
independent of the relative positions and orientations of the objects. 
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C. Euclidean Electromagnetic Action 

By replacing the time T by —ihjS, we transform the Minkowski space functional integral 
Z{T) into the partition function Z{P). In = gauge, the result is simply to replace the 



frequencies Un = ^ in Eq. (II. 4) by = icK^, where /€„ is the n^^ Matsubara frequency 
divided by c. (In other gauges the temporal component of the vector field must be 
rotated too.) 

The Lagrangian is quadratic, so the modes with different /€„ decouple and the partition 
function decomposes into a product of partition functions for each mode. Since the elec- 
tromagnetic field is real, we have Fj*{uj) = E(— cj) on the real axis. We can thus further 
simplify this decomposition on the imaginary axis by considering k, > only, but allowing 
E and E* to vary independently in the path integral. Restricting to positive k is possible 
because the response functions e{icK,x.) and fi{icK,,x) are invariant under a change of sign 
in icK, as shown Appendix |A| In the limit P ^ oo, the sum ^„>o turns into an integral 
^ Jq°° dn, and we have 

he 

So = -— / dn logZ(ft), (11.11) 



where 

Z{k) = I VAVA* exp 



27r 



-p I dxE* ■ ( 1+ X Vx ) E + ^E* ■ V(zc/t,x)E 



(11.12) 

V(2CK,x) = I/t2(e(icK,x) - 1) + V X ( — ^ — - - 1 ) V x . (11.13) 

\/i(zcfi:, x) J 

The potential Y{icK, x) is real for real k, even though e and fi can have imaginary parts for 
real frequencies uj. Our goal is now to manipulate Z[k) in Eq. (11.12) so that it is computable 
from the scattering properties of the objects. 

III. GREEN'S FUNCTION EXPANSIONS AND TRANSLATION FORMULAS 

A. The free Green's function 

The free Green's function and its representations in various coordinate systems are crucial 
to our formalism. The free electromagnetic field (V = 0) obeys equations of motion obtained 



by extremizing the corresponding action, Eq. (II. 2) 

^2 



+ V X Vx ) E(cj,x) = 0. (III.l) 



We will employ the electromagnetic dyadic Green's function Go, defined by 

(^-I^ + V xVx^Goiiu, X, x') = M^^) (x - X ) , (III.2) 

written here in the position space representation. It is easy to express Gq as a Fourier 
transform, 

r fllr pik-{x-x') / \ 
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where the displacement of the singularities at /c = ±^ corresponds to outgoing wave bound- 
ary conditions at infinity. By replacing the factors of k by gradients, Gq may be expressed 
in terms of elementary functions, 



.o(cc;,x,x 



tu^ / 47r|x 



^iu!\x—x.'\/c 



(III.4) 



In this representation it is easy to see that Go is transverse, i. e. V-Go(x, x', u) = 
V = 0, for X 7^ x'. Go is not transverse at x = x', as can be seen by taking the divergence 



of Eq. (III.2). 



We work in coordinate systems in which we can use separation of variables and employ a 
spectral representation of Go(x, x',a;). That is, we represent the Green's function through 
the complete set of non-singular ("regular"), transverse solutions to the differential equation. 



Eq. (III.l) 



E 



reg. 



x) = (x|E-^(a;)), 



(III.5) 



represented formally by the eigenstate kets \ W^^{u))), where the generalized index a labels the 
scattering channel, including the polarization. For example, for spherical wave functions it 
represents the angular momentum quantum numbers (/, m) and the polarization E or M. We 
will choose to normalize these states in accord with standard conventions in electromagnetic 
scattering theory; as a result they are not necessarily normalized according to the conventions 
typically used in quantum mechanics. A list of the eigenfunctions for various common bases 
is given in Appendix |Bj The Green's functions can be expressed as the coordinate-space 
matrix element of the operator 



/■oo 

Jo 



[UJ 



(III.6) 



where the ie has again been included to implement outgoing wave boundary conditions, so 
that the Green's function is causal.^ We use the symbol Gq to represent both the matrix- 



valued representation of the Green's function in position space, Eq. (III.2), and the abstract 



Hilbert space operator, Eq. (III. 6). The coefficients Ca{uj') are inserted because of our choice 



of normalization and ensure that 



POO 

JO 



)I=I- 



(111.7) 



It is also useful to represent the Green's function in a different way, in which one of 
the separable coordinates is identified as the "radial" variable and treated differently from 
the remaining coordinates. We let represent this coordinate and denote the remaining 
coordinates as ^2 and ^3. We introduce the "outgoing" solution in ^1, which is in the same 



The coordinate space matrix element of Eq. (III. 6) is transverse for all x and x', and therefore differs 



from the Green's function defined in Eq. (III. 4 1 by terms local at x = x'. Since we never employ Gq at 



coincident points, we ignore this subtlety |73j . The use of the retarded Green's function not only makes 
sense physically, but is also dictated by the imaginary- frequency formalism, just as is the case for the 
response functions e and /i. It is the retarded response functions that are analytically continued in the 
frequency domain to positive imaginary frequency, as shown in Appendix 
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scattering channel as the corresponding regular solution but obeys outgoing wave boundary 
conditions as — >■ oo. It is linearly independent of the regular solution. The full outgoing 
solution is then obtained by multiplying the outgoing solution for by the regular solutions 
for ^2 and ^3. We can then express one of the regular wave functions in the position space 



representation of Eq. (III. 6 ) as a sum of the outgoing solution for uj and the outgoing solution 
for —uj. By specifying explicitly which of the two arguments of the Green's function has a 
greater value of .^1, we can carry out the uj integral for each of these two terms separately 
by closing the contour in the appropriate half-plane [73] , and obtain 

In this form, the outgoing wave boundary condition is implemented explicitly. Since the 
Green's function is written as a linear combination of solutions to the free wave equation, 



it clearly satisfies Eq. (III. 2) for x 7^ x'. The normalization constant Ca{uj), which is 
determined using the Wronskian of the regular and outgoing solutions and the completeness 
relationship for the regular solutions in ^2 and ^3, sets the correct "jump condition" for 



X = x'. 



The outgoing solution is typically singular at ^1 = 0, but the Green's function with 
distinct arguments does not encounter that region, because the outgoing function is always 
evaluated for the larger argument. For example, in a spherical system the outgoing solution 
could take the form of a spherical Hankel function h^^\kr) ~ with k = 00 /c, which 
obeys outgoing wave boundary conditions, is singular at the origin, and is independent of 
the corresponding regular solution ji{kr). 

We will usually work on the imaginary k-axis, in which case we will encounter the cor- 
responding modified special functions. We continue to label these functions as "regular," 
"outgoing," and "incoming," even though they now increase exponentially for large ^1 for 
incoming and regular waves and decrease exponentially for outgoing waves. We also note 
that it may be convenient to redefine the wave functions to match the standard form of the 
corresponding modified functions, and to assign different phases to the two polarizations. 
The prefactor Ca{w) is then correspondingly redefined as Cai^) to incorporate these changes. 
A list of Green's function expansions in various common bases is given in Appendix [B] 

For a Cartesian coordinate system some of the previous statements have to be adapted 
slightly. We will take one of the Cartesian coordinates, say z, to be the "radial" coordinate, as 
required by the context. For example, z might be the direction normal to the planar surface 
of a dielectric. The solutions are then given in terms of plane waves, Q^'^^^+^'^yV^^V^'^/'^'i^'^^^ ^ 
where is the momentum perpendicular to the z direction. All are regular and all con- 



tribute in the integral representation of Eq. (III.6). After analytic continuation to imaginary 



frequency, the free Green's function in Cartesian coordinates is expressed by the above for- 
mula if we identify outgoing solutions with plane wave functions that are exponentially 
decreasing in the +z direction, q^^^^+'^^vv- V'^^+^'i ^ ^ and regular solutions with the exponen- 
tially growing solutions e^^^-^j+ifc^y+^/i^+i^^^ 

The wave functions that appear in the series expansion of the free Green's functions 



in Eq. (III. 8) satisfy wave equations with frequency uo. The integral representations in 



Eq. (III.6), on the other hand, contain wave functions of all frequencies. As we will see in 



Sect. IV, the ability to express the Casimir energy entirely in terms of an "on-shell" partial 



wave expansion with fixed u will greatly simplify our calculations. 
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B. Translation matrices 



We will use the free Green's function described in the previous subsection to combine the 
scattering amplitudes for two different objects. In this calculation, the one argument of the 



Green's function will be located on each object. As a result, if Eq. (III.8) is written in the 
basis appropriate to one object, we will want to "translate" one of the scattering solutions 
to the basis appropriate to the other object. The configuration of the two objects — either 
outside of each other, or one inside the other — will determine which object has the larger 
or smaller value of ^i, and therefore which solution needs to be expanded in the other basis. 
We will make use of two expansions: 

1. The regular solutions form a complete set no matter what origin is used to define the 
decomposition into partial waves. Let {E|J3^^(k, Xj)} be the regular solutions expressed 
with respect to the origin of coordinates appropriate to object j, Oj. It must be 
possible to expand a regular solution E^''^(K,Xj), defined with respect to the origin Oi 
appropriate to object z, in terms of the {E|J3''^(k, Xj)}, 

E-^(«:,x,) = ^ V^;j«:,X,,)E-^(«:,x,), (III.9) 

a 

where Xjj = — Xjj = Xj — x^ is shown in Fig. [Tj Note that Xj and Xj refer to the same 
space point x, expressed as the displacement from different origins. This expansion 
will be applicable to the case of one object inside the other. 

2. The same type of expansion must also exist when the original wave obeys outgoing 
boundary conditions except in a region that contains the origin Oi, where E°"*(/t, Xj) 
is singular. We therefore have the expansion 

Er(«:,xO = ^W^;Jk,X,OE-^(«:,x,), forx ^ N{Oi) (III.IO) 

where N{Oi) is a neighborhood of the origin Oi. This expansion will be applicable to 
the case where the objects are outside each other. 

To apply these results to a given geometry, we must be able to distinguish between regular 
and outgoing waves over the whole of each object. That is, we require there to exist an origin 
and a separable coordinate system so that for all points x in one object and x' in another 
object, ^i(x) is always greater than ^i(x'), or vice versa. Having ^i(x) > ^i(x') ensures 
that the Green's function is always evaluated by letting x be the argument of the outgoing 
wave function and x' be the argument of the regular wave function. We therefore require 
that any two objects be separated by a surface defined by the locations x where ^i(x) is 
constant, as shown in Fig. [1} Depending on the coordinate system, this surface could be a 
plane, cyhnder, sphere, etc. 

The case of an elliptic cylinder and a circular cylinder illustrates this requirement. At 
large distances, the elliptic cylinder object can be separated from the circular cylinder object 
by a circular cylinder of radius p, as shown in Fig. [2^. All points on the elliptic cylinder 
object have values of pi that are smaller than any point on the circular cylinder object, so 
in this case we could carry out the calculation in ordinary cylindrical coordinates. However, 
as shown in Fig. [2]d, if the separation becomes small enough, points on the circular cylinder 
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0\ 





FIG. 1: (Color online) Geometry of the outside (left) and inside (right) configurations. The dotted 
lines show surfaces separating the objects on which the radial variable is constant. The translation 
vector Xjj = Xj — Xj = — X^j describes the relative positions of the two origins. 



object are closer to the center of the elliptic cylinder object {i.e. they lie at smaller pi than 
points on the elliptic cylinder object), and our method cannot be used in ordinary cylindrical 
coordinates. However, in elliptic cylindrical coordinates (see Appendix B4), the surface of 
the elliptic cylinder object is itself a surface of constant elliptical radius /ii, so all points 
on the elliptic cylinder object have smaller /ii than any point on the the circular cylinder 
object, and our method applies. This case is shown in Fig. 

In a plane wave basis, we would exclude the case of two interlocking combs [71], since 
each comb has values of z that are both bigger and smaller than points on the other object, 
so again a single assignment of regular and outgoing solutions cannot be made. 

When object j lies wholly outside of object i, as shown in the left panel of Fig. [1} in the 
basis of object i the point on object j will always have greater than the point on object 
i. We will therefore need to expand the outgoing wave in the basis for object j. Since the 
origin Oi is never encountered when the point x lies on object j, the outgoing solutions for 
i can be expanded in terms of the regular solutions for object j using Eq. (III. 10). Since i is 
also wholly outside j, we can also proceed the other way around and expand the outgoing 
wave functions in the basis of object j in terms of regular solutions in the basis of object 
i. This implies that the translation matrix satisfies = U^'^''^ . When one object is inside 
another, as shown in the right panel of Fig. [T| in the basis of object z, the point on object 
j will always have smaller than the point on object i. We will therefore need to expand 
the regular wave in the basis for object j using Eq. (III.9). In contrast, we cannot use the 
expansion of the outgoing wave functions, because the origin of the inside object may overlap 
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( (a) (b) (c) 



FIG. 2: An elliptic cylinder approaching another cylinder. When the elliptic cylinder is far (a), 
every point on the cylinder has smaller radius than any point on the lower cylinder and an expansion 
using an ordinary cylindrical basis can be used. This expansion breaks down once the elliptic 
cylinder is close (b), but in that case an expansion using an elliptic cylindrical basis applies (c). 



with the origin of the outside object, in which case the expansion does not converge. 

For a Cartesian geometry, the translation matrix is proportional to 
g-«kx-x^i,x- Vt^+k^^ii,z _ takes this simple form because plane wave functions are 
eigenfunctions of the translation operator. Then the "regular" wave function is evaluated 
on the object whose z coordinates are smaller and the outer and inner objects have larger 
and smaller z values, respectively. 

The criterion for the expansion of the outgoing or regular wave functions is not topological. 
Instead, the proximity of the objects and their origins determines which expansion to use. 
In practice, it is usually easy to see which expansion is appropriate for any objects. 

After expanding wave functions with respect to another origin using translation matrices, 
we can convert the wave functions from one basis to another, for example from plane wave to 
spherical or cylindrical wave functions. This transformation is useful when the two objects 
are best described in different coordinate bases. The needed conversion matrices are supplied 
in Appendix [Dj Since it is more convenient to describe this conversion as a change of basis 
of the scattering amplitudes, we will not explicitly consider the combination of translation 
and conversion in this derivation, but instead we will illustrate the change of basis of the 
scattering amplitude in the examples. 

C. Green's functions and translation matrices 

To obtain the Green's function when one argument, say x, lies on object i and the other 
argument, say x', lies on object j, we expand Go^icn, x, x') in terms of coordinates Xj and x^- 
that describe each point relative to the origin of the body on which it lies. For the different 
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situations given above we have 
GolicK, X, x') = 

EJJfs(fi;,Xi) W;^^(K)E^''^*(fi:, Xj) if i and j are outside each other 



E^^s(fi:, Xi) ® V^p{K)¥.'^* (k, x;.) I if ^ ^gi^^ J ' (plane wave basis) 

Er(«:,xO ® w^Uk)e^^*^*(k,x;.) I j^^'f'^'^ ':^\ , . . 

" ^ ' apy / ij ^ ' J ^ It J IS below « (plane wave basis) 



(III.ll) 



where W;^*^(k) = V^^^ {k) and Co is the normalization constant defined in Eq. (III.8). 
We can express these cases in the consolidated form, 



Go(zc«:,x,x') = Y.C,{n) (E-^(«:,x.) Er(«:,x.)) ® [^^f^l^ ^^^i^^ ) 

(III.12) 

where only one of the three submatrices is nonzero for any pair of objects i and j as given 



in Eq. (III.ll). The expansion can be written more formally as 



a,f3 



(Ei?(«:) 



(III.13) 



where the bras and kets are to be evaluated in position space in the appropriately restricted 
domains and the X matrix is defined, for convenience, as the negative of the matrix contain- 
ing the translation matrices. 

The translation matrices for various geometries are provided in Appendix O 



IV. A REVIEW OF ASPECTS OF THE CLASSICAL SCATTERING OF 

ELECTROMAGNETIC FIELDS 

In this section, we review the key results from scattering theory needed to compute the 
scattering amplitude of each body individually. Comprehensive derivations can be found in 
Refs. [751 EB] ■ The approach we will use was first developed by Waterman [771 CH] , albeit 
not in the operator form that is used here. In the subsequent section we will then combine 
these results with the translation matrices of the previous section to compute Z{k). 

A. Electromagnetic scattering 



The Fourier-transformed electromagnetic action of Eq. (II. 2) yields the frequency- 
dependent Maxwell equations: 



V X E(co',x) = i—B(uj,x.) 

c 



V X -B(cj,x) 



UJ 

-i—eE[u!, x). 
c 



(IV. 1) 
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Combining these two equations, we obtain 



where 



V(cu, x))E(a;, x) = ^^E(u;, x), (IV.2) 
Ho = V X Vx, 

(IV.3) 



V(.,x)^I-(l-e(.,x)) + Vx^— -iJVx, 

which is the same potential operator as the one obtained by rearranging the Lagrangian (see 
Eq. (II.4)). Since the electromagnetic potential is a differential operator, care must be taken 
with operator ordering. 

The Lippmann-Schwinger equation [79] 

|E) = |Eo) - GoV|E) (IV.4) 



expresses the general solution to Eq. (IV.2). Here Gq is the free electromagnetic ten- 



sor Green's function discussed in Sec. Ill and the homogeneous solution |Eo) obeys 

^— ^I + Hoj |Eo) = 0, which can be chosen to be either a regular or outgoing wave for 

a particular frequency ui. We can iteratively substitute for |E) in Eq. (IV.4) to obtain the 
formal expansion 

|E) = |Eo) - GoV|Eo) + GoVGoV|E) - . . . 
= |Eo) — GoT|Eo), 
where the electromagnetic T-operator is defined as 

and G is the Green's function of the full Hamiltonian, ( — + Hn + V ) G = I. We note 



that T, Go, and G are all functions of frequency uj and non-local in space. As can be seen 
from expanding T in Eq. (IV. 6) in a power series, T(a;,x, x') = (x|T(a;)|x') is zero whenever 
X or x' are not located on an object, i.e., where V(u;,x) is zero. This result does not, 
however, apply to 

T-i = Go+V"\ (IV. 7) 

because the free Green's function is nonlocal. 

Next we connect the matrix elements of the T-operator between states with equal uo to the 
scattering amplitude JF. In our formalism, only this restricted subset of T-operator matrix 
elements is needed in the computation of the Casimir energy. 



1. Outside scattering 

We consider a scattering process in which a regular wave interacts with an object and 
scatters outward, as depicted in the left panel of Fig. |3}^ For outside scattering the homo- 



Alternatively, we can set up asymptotically incoming and outgoing waves on the outside and regular 
waves inside. The amplitudes of the outgoing waves are then given by the iS-matrix, which is related to 
the scattering amplitude T hy T = {S — X)/2. Although these two matrices carry equivalent information, 
the scattering amplitude will be more convenient for our calculation. 
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FIG. 3: (Color online) The scattering waves for outside scattering (left panel) and inside scatter- 
ing (right panel). In both cases the homogeneous solution Eo(cj) is shown in bold. For outside 
scattering, the homogeneous solution is a regular wave, which produces a regular wave inside the 
object and an outgoing wave outside the object. For inside scattering, the homogeneous solution is 
an outgoing wave, which produces a regular wave inside the object and an outgoing wave outside 
the object. 



geneous solution |Eo) in Eq. (IV.5) is taken to be the regular wave function \W^^{uj)). We 
choose a convenient "scattering origin" in the inside region, consistent with any symmetries 
of the problem if possible. 

To find the field E at a coordinate x far enough outside the object, we use Eq. (IV.5) in 
position space and the expansion in Eq. (III. 8) for Gq: 



E(w,x) 



E 



UJ, X 



- J2 Er(^' ^) / C(s{uj)E'^^*{uj, x') ■ T(cu, x', ^")E'^^{uj, x")rfx'rfx". 

(IV.8) 

Here "far enough outside" means that x has larger than any point on the object, meaning 
that we are always taking the same choice in Eq. (III.8), as described in Sect. Ill The 
equation can be written in Dirac notation, again with the condition that the domain of the 
functional Hilbert space is chosen appropriately to the type of solution. 



|E(^)) 



)) + Y.\'^Ti^)) X (-l)C,(a;)(E-^(c.)|T(a;)|E-^(a;)), 

(3 ' 



(IV.9) 



which defines as the exterior /exterior scattering amplitude (the one evaluated between 
two regular solutions). We will use analogous notation in the other cases below. 

At coordinates x "far enough inside" a cavity of the object, meaning that x has smaller 
,^1 than any point on the object, we have the opposite case in Eq. (III. 8) and the field E is 
given by 



\E{u:)) = \E-^{u)) + J2 lE^M) X (-l)C7,(a;)(E^^(a;)|T(c.)|E-^(c.)), 

(3 ' 



(IV.IO) 



where again the free states are only defined over the appropriate domain in position space, 
and JF*^ indicates the interior /exterior scattering amplitude. 



16 



2. Inside scattering 

In the study of Casimir problems with one object inside the other, it is useful to imagine 
a situation that would be difficult to realize in actual scattering experiments, in which the 
wave probing the object originates inside the object and is scattered as a regular wave inside 
the object and as an outgoing wave outside, as depicted in the right panel of Fig. |3| 

The situation is expressed mathematically by letting the homogeneous solution |Eo) in 



Eq. (IV. 5) be an outgoing wave |E°'^*(ci;)). The equation can be expressed in condensed form 



as before. Inside the object we have 

|EM) = |Er H) + J2 lE^M) X (-i)c^/3M(E;?M|TH|Er H), (iv.ii) 

and outside the object we have 

|EH) = |Er M) + l^rH) X (-l)C,H(E-^H|TM|Er(c.)) . (IV.12) 

/3 ^ ' 



3. Remarks 

We have obtained the scattering amplitude in the basis of free solutions with fixed u. 
Since one is normally interested in the scattering of waves outside the object, the scattering 
amplitude usually refers to JF*^^. We will use a more general definition, which encompasses all 
possible combinations of inside and outside. The scattering amplitude is always "on-shell," 
because the frequencies of the bra and ket wave functions are both equal to u. As a result, 
it is a special case of the T-operator, which connects wave functions with different u. 

It is usually not practical to calculate the matrix elements by finding the abstract T- 
operator and taking its inner products with free wave functions. Instead, one typically 
considers an ansatz for the solutions appropriate for inside or outside scattering in the 
various regions, with unknown scattering amplitudes, and then solves the wave equation, 
matching the solutions in different regions at their boundaries. 

We will find it convenient to assemble the scattering amplitudes for inside and outside 
into a single matrix. 



{-1)C^{k) 



{W^^{K)\T{tcK)\Y.';\K)) (E-s(«:)|T(^c«:)|E-*(«:)) 
(E-(«:)|T(zm)|E^^^(«:)) {Y}-{^)\T{ic^)\Ef\^)) 



(IV.13) 



Here we have written this expression in terms of modified wave functions for uo = icK, 
with the corresponding normalization constant, since that is the case we will use. The 
derivations of the scattering amplitudes carry over directly to this case, with k replaced by 
uj; for example, Eq. ( IV.9[ ) becomes 



(IV.14) 
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V. PARTITION FUNCTION IN TERMS OF THE SCATTERING AMPLITUDE 



With the tools of the previous two sections, we are now able to re-express the Euclidean 
electromagnetic partition function of Eq. (11.12) in terms of the scattering theory results 
derived in Section IV for imaginary frequency. We will exchange the fluctuating field A, 
which is subject to the potential V(2ck, x), for a free field A', together with fluctuating 
currents J and charges — ■ J that are confined to the objects. The sequence of two 
changes of variables that will be performed is often referred to as the Hubbard-Stratonovich 
transformation in condensed matter physics. 



We multiply and divide the partition function Eq. (|II.12|) by 



j rfxJ*(x) • V-^(icfi:,x)J(x) 



detV(icK,x,x') , (V.l) 



where l^^^j indicates that the currents are defined only over the objects, i.e. the domain 
where V is nonzero (and therefore exists), and we have represented the local potential 
as a matrix in position space, Y{icK, x, x') = Y{icK, x)5*^^^(x— x'). Our derivation generalizes 
straightforwardly to the case of a nonlocal potential V(ic/t, x, x'), assuming it is still confined 
to each object individually. 

We then change variables in the integration, J(x) = J'(x) + -V(icK, x)E(x) and J*(x) = 
j'*(x) + iV(icK,x)E*(x), to obtain 



Z(k) 



1 

W 



VAVA* VJ'VJ' 



I obj 



exp 



T 
H 



rfx n 



K 



+ J'*(x) + -V(zcK,x)E*(x) • V"X^c«,x) J'(x) + -V(icft,x)E(x 



K 



(V.2) 



where 



7^ = E* (x) ■ I + — V X V X E(x) + — E* (x) • N{icK, x)E(x) 



(V.3) 



Next we use a second change of variables, E(icK, x) = E'(icK, x) — m J (hd GQ^icn, x, x^) J^(x^) 



W 



VJ'VJ' 



I obj 



and E*{icK,x.) = E'*(icK,x) — in J fix' Go(«ck, x, x') J'*(x'), which simplifies Eq. (V.2) to 

7. 

Z{k) 

exp 

where 



T 

li 



j c/xrfx' J'*(x) • (Go(icK,x,x') + V-^(ic/t,x,x')) J'(x') 



(V.4) 



VAIVA'* exp 



T 

li 



dxE'*(x) • 1+ X Vx E'(x) 



(V.5) 



is the partition function of the free field, which is independent of the objects. The new 
partition function of Eq. (V.4) contains a sum over current fluctuations in place of the 
original field fluctuations in Eq. (11.12). The interaction of current fluctuations at different 
points X and x' is described by the free Green's function Go(icK, x, x') alone. (If the potential 
Y{icK, X, x') is nonlocal, this statement still holds for two points x and x' on two different 
objects.) This is the expected interaction term. For example, in the static limit k = 0, the 
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free Green's function is just the Coulomb interaction term ^^|^^_^,| . The inverse potential 
penalizes current fluctuations if the potential is small. In vacuum, the potential vanishes, 
so current fluctuations are infinitely costly and thus are not permitted. But of course the 
current fluctuations are already constrained to the objects. 

To put the partition function into a suitable form for practical computations, we will 
use the results of the previous sections to re-express the microscopic current fluctuations as 
macroscopic multipole fluctuations, which then can be connected to the individual objects' 
scattering amplitudes. This transformation comes about naturally once the current fluctu- 
ations are decomposed according to the objects on which they occur and the appropriate 
expansions of the Green's function are introduced. We begin this process by noticing that 



the operator in the exponent of the integrand in Eq. (V.4) is the negative of the inverse of 



the T-operator (see Eq. (IV.7)), and hence 



Zn detV 



~^(iCK, X, x'l 



det T{icK, x, x') 



(V.6) 



which is in agreement with a more direct calculation: Since Zq = det Go(^c/t, x, x') and 
Z{k) = det G(zcfi;, X, x'), we only need to take the determinant of Eq. (IV. 6) to arrive at the 
result of Eq. (V.6). 

Both Zq and det V~^(icK, x) are independent of the separation of the objects, since the 
former is simply the free Green's function, while the latter is diagonal in x. Even a nonlocal 
potential V(ic/t, x, x') only connects points within the same object, so its determinant is 
also independent of the objects' separation. Because these determinants do not depend on 
separation, they will be canceled by a reference partition function in the final result. We are 
thus left with the task of computing the determinant of the T-operator. 



A. From the determinant of the T-operator to the scattering amplitude 



As has been discussed in Sec. IV, the T-operator T(2ck, x, x') is not diagonal in the 
spatial coordinates. Its determinant needs to be taken over the spatial indices x and x', 
which are restricted to the objects because the fluctuating currents J(x) in the functional 
integrals are zero away from the objects. This determinant also runs over the ordinary vector 
components of the electromagnetic T operator. 

A change of basis to momentum space does not help in computing the determinant of the 
T-operator, even though it does help in finding the determinant of the free Green's function. 
One reason is that the momentum basis is not orthogonal over the domain of the indices x 
and x', which is restricted to the objects. In addition, a complete momentum basis includes 
not only all directions of the momentum vector, but also all magnitudes of the momenta. So, 
in the matrix element (Ek|T(ci;)|Ek') the wave numbers k and k' would not have to match, 
and could also differ from uj/c. That is, the matrix elements could be "off-shell." Therefore, 
the T-operator could not simply be treated as if it was the scattering amplitude, which is 



the on-shell representation of the operator in the subbasis of frequency u (see Sec. IV), and 
is significantly easier to calculate. Nonetheless, we will see that it is possible to express the 
Casimir energy in terms of the on-shell operator only, by remaining in the position basis. 



From Eq. (IV.6), we know that the inverse of the T-operator equals the sum of the free 
Green's function and the inverse of the potential. Since the determinant of the inverse 
operator is the reciprocal of the determinant, it is expedient to start with the inverse T- 
operator. We then separate the basis involving all the objects into blocks for the n objects. 
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In a schematic notation, we have 



|T-^|x')] 



[(xi|Tr^ 


Ix'i)] 


[(xilGol 


x'2)] 




^(x2|Go 


X'l)] 


[(X2|T2^ 


lx'2)] 











(V.7) 



where the ij^^ submatrix refers to x G object i and x' G object j and Xj represents a point 
in object i measured with respect to some fixed coordinate system. Unhke the position 



vectors in Sec. |III[ at this point the subscript of Xj does not indicate the origin with respect 
to which the vector is measured, but rather the object on which the point hes. Square 
brackets are used to remind us that we are considering the entire matrix or submatrix and 
not a single matrix element. We note that the operators T and Go are functions of icn, but 
for simplicity we suppress this argument throughout this derivation. When the two spatial 
indices lie on different objects, only the free Green's function remains in the off-diagonal 
submatrices. Even if the potential Y{icK, x, x') is nonlocal in space, it does not connect 
points on different objects. It follows that the inverse of the potential is block diagonal in 
position space, where each block involves points on the same object, i.e., (xj|V~^|x^) = 
for i j. 

Next, we multiply by a reference T-operator Too without off-diagonal submatrices, 
which can be interpreted as the T-operator at infinite separation. 



[(x|TooT-^|x")] 



[(xi 


|x'/)] 


[/ rfx;(xi|Ti| 


x;)(x;|Go|x'2')] 




[/ rfx'2 (X2IT2I 


x'2)(x'2|Go|x;')] 


[(X2 













(V.8) 



Each off-diagonal submatrix [J (ix^(xj|Tj|x^)(x'j|Go|xJ)] is the product of the T-operator of 
object i, evaluated at two points Xj and x'^ on that object, multiplied by the free Green's 
function, which connects x'^ to some point x" on object j. 

Now we shift all variables to the coordinate systems of the objects on which they lie. As a 
result, the index on a position vector Xj now refers to the object i on which the point lies and 
to the coordinate system with origin Oi in which the vector is represented, in agreement with 
the notation of Sec. Ill The off-diagonal submatrices in Eq. (V.8) can then be rewritten 



using Eq. (III. 13) as. 



E 

a,l3 



((x,|T,|E-^(«:)) (x,|T,|Er(/.)))X:j,| ^^^^^ ] {-C,{^)) 



(V.9) 



The matrix [(x|TooT |x")] has the structure I + AB, where 



a 






[(xi|Ti|E-g(«:))Xa [(xi|Ti|Er(«:))Xia 




[(x2|T2|E-s(«:))X2y [(x2|T2|Er('^))X^y 














(V.IO) 
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and 



( [-C,{k){Y.-\k)\^{)] 
\-C,{k){Y}^{k)\^{)] 








V 








rcg. 



K X 



K x: 



)] 
^2)] 



(V.ll) 



and the matrix multiplication now encompasses both the object index and the partial wave 
index (3. Although the same symbols are used for each wave function, the bases (spherical, 
planar, etc.) can be chosen differently for each object. 

Using Sylvester's determinant formula det(I + AB) = det(I + BA), we see that the deter- 



minant is unchanged if we replace the off-diagonal submatrices in Eq. (V.8) by 



|T,|E-*(«:)) 



(V.12) 



With this change, the diagonal submatrices in Eq. (V.8) become diagonal in the partial 
wave indices rather than in position space. The matrix elements of the T-operator are 
the scattering amplitudes, which can be obtained from ordinary scattering calculations, as 
demonstrated in Sec. IV The first matrix in Eq. (V.12), including the prefactor (— 1)Cq,(k), 



is Fj(K), the modified scattering amplitude of object i, defined in Eq. (IV. 13) 



Putting together Eqs. (11.11), (11.12), (V.6), and (V.8), we obtain 



he 
2^ 



where 



dfi:logdet(MM^^), 



Fr' 

X21 X23 



(V.13) 



(V.14) 



and Mq^^ is a block diagonal matrix diag(Fi F2 
Using the block determinant identity 



det 



C 



det (A) det (D - CA^^B) = det (D) det (A - BD^^C) , (V.15) 



we can simplify this expression for the case of the interaction between two objects. 



he 
2^ 



d/tlogdet (I - F^X^'^'FtX'"') 



(V.16) 



Usually, not all of the submatrices of F and X are actually needed for a computation. For 
example, if all objects are outside of one another, only the submatrices T*^^ of the scattering 
amplitude that describe outside reflection are needed. If there are only two objects, one 
inside another, then only the inside reflection submatrix JF** of the outside object and the 
outside reflection submatrix T^'^ of the inside object are needed. 

In order to obtain the free energy at nonzero temperature instead of the ground state 
energy, we do not take the limit /? ^ 00 in Eq. (II. 9) [23J. Instead, the integral ^ dtv 
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is replaced everywhere by ^ where ck„ = ^ with n = 0, 1, 2, 3 ... is the nth Matsub- 
ara frequency. A careful analysis of the derivation shows that the zero frequency mode is 
weighted by 1/2 compared to the rest of the terms in the sum; this modification of the sum 
is denoted by a prime on the summation symbol. The factor of 1/2 comes about because 
the fluctuating charges or currents have to be real for zero frequency. Thus, for kq, the 
expressions on the right hand side of Eq. (V.6) should be placed under a square root. (For 
a complex field, both signs of the integer n would be included separately, and n = would 
be included once, with the normal weight.) 

If the medium between the objects is not vacuum but instead has permittivity tm{icK) 
and magnetic permeability /im(^CK) different from unity, then the free Green's function 
is multiplied by ^miicn), and its argument n is replaced by where 
■\/ em{icK) i^rniicn) is the medium's index of refraction. Effectively, this change just scales 
all frequency dependencies in the translation matrices X(k), which become li{nm{icK)K). 
Furthermore, the scattering amplitudes absorb the factor fim{icK,) from the free Green's 
function and change non-trivially, i.e. not just by some overall factor or a scaling of the 
frequency. They have to be computed with the nonzero electric and magnetic susceptibilities 
of the medium. 



VI. APPLICATIONS 

Our technique for calculating the Casimir energy applies to a wide range of situations. 
In this section we demonstrate the method through a variety of examples. 



A. London and Casimir-Polder interaction between two atoms 

As a simple example, we re-derive the interaction between two identical neutral atoms in 
the ground state [19]. The atoms are described in a two-state approximation. Within this 
approximation, the electric dipole polarizability of the atoms is given by 



2 f 

a = 2 ^' 

where e is the electron charge, m is the mass, /lo is the oscillator strength of the 0^1 
transition, and ujiq is the frequency of that transition. We perform a Wick rotation uj — *■ 
icK and set k, = u/d, where d is the distance between the atoms. By introducing the 
characteristic length scale dio = c/uio and the static electric polarizability 

"0 = fiorodlo , (VI.2) 

where vq = e^/(mc^) ^ 10~^^m is the classical electron radius, the polarizability can be 
written as 

{d/dio)'^ao 
{d/dioY + 

In the isotropic-dipole approximation, the only nonzero element of the scattering ampli- 
tude of the atom is given in terms of the electric dipole polarizability as 



22 



for m = —1, 0, 1. The atoms are assumed to have no magnetic polarizabihty. Using the 



general expression of Eq. (V.16) for the interaction energy between two objects, we get 



S 



he 

27rd 



du log 



l-4(l + n)V2"-^ 




1 - (1 + « + V^'" 



■-^ f du {a'^{u)Y (3 + 6u + hu^ + 2u^ + u^) e 
nd Jq 



2u 



(VI.5) 
(VI.6) 



where we have expanded the logarithm assuming a^{u) <C d'^, so that the interaction energy 
is proportional to the product of the polarizabilities of the atoms. It is instructive to consider 
two limits. First, assume that d <^ diQ. By the change of variable u = {d/dio)z, one easily 
finds that in Eq. (VI.6), only the leading constant term of the polynomial in u has to be 



considered, and the exponential factor can be ignored. The integral is convergent at large 
u due to the behavior of the polarizabihty a^{u) at large u. The integral over z yields, to 
leading order in d/dio, the energy 



St 



3, «o 
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which is the well-known London interaction 

In the opposite limit d ^ rfio retardation is important. From Eq. (VI. 3) we see that the 
frequency (m) dependence of the polarizabihty now can be neglected, so that a ~ oq. In 
this retarded limit, the polynomial and exponential in u in Eq. (VI.6) are important, and 
integration yields the energy 

P 23 a2 

which is known as the Casimir-Polder interaction 



(VI.8) 



For general distances d, the interaction between the atoms can be computed numerically; 
to quadratic order in the polarizabihty it is given by the integral of Eq. ( VI.6 ). The numerical 



result and the two limiting forms of the interaction are shown in Fig. 4[ 



B. Derivation of the Lifshitz formula 



Next we consider two semi-infinite half-spaces with uniform electric and magnetic sus- 
ceptibility, as depicted in Fig. [s] |23l |2ll |25l |26l [27]. We choose a plane wave basis oriented 
along the z axis. 

We decompose the scattering amplitude into magnetic (transverse electric) modes M and 
electric (transverse magnetic) modes N. 

For the upper object a the scattering solution is [81] 



E(«:,x) = N-*(/s:,x) + 



L'^dk'^ 
(27r)2 

L^d^k'^ 
(27r)2 



M,^,^(k, x)JF**,^, jv/,kxAf + -l^k'^l^' ^)-^aV| E,l<ij_ 



M 



(VI.9) 



Here L is the length of each side of the plates, is the momentum perpendicular to the z 
direction, and the subscripts M and E on the scattering amplitudes denote the magnetic and 
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FIG. 4: (Color online) Interaction energy £ of two identical atoms, Eq. (VI.6), as a function of their 
separation d. The curve shows the crossover between the London interaction {d <C diQ = c/ljio), 



Eq. (VI. 7), and the Casimir-Polder interaction (d 3> dio); Eq. (VI. 8). 




FIG. 5: (Color online) The upper infinite half space a is located a distance d above the half space 
h. This is the original configuration considered by Lifshitz. Each half space has its own uniform 
electric permittivity ei{icK) and magnetic permeability ^i{icK). We note that our calculation holds 
even if the two origins Oa and Oh are displaced horizontally from one another, as shown here. 



electric polarizations respectively. We consider the limit L ^ oo. The scattering amplitudes 
are given by 
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in terms of the Fresnel coefficients 
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Here, is the index of refraction, na{icK) = ^Jea{icK)^Xa{'i'CK). In the hterature the Fresnel 
coefficients are also sometimes labeled with s instead of M and p in place of E. 

The lower object h has the same scattering properties. The relevant scattering equation 
is the same as in Eq. (VI. 9), with "reg" and "out" exchanged and JF" replaced by T^^, 



which is obtained from JF** simply by substituting the permittivity eb{icK) and permeability 
fib{icK) in place of those of object a. 

Using the appropriate X submatrices as specified in Eq. ( III.11[ ) and the corresponding 
submatrices of F, the energy (V.16) for two objects can be expressed as 



he 
2^ 



POO 

/ dK log det (T - jr^^iy\;bajree^ba^ 
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Since the matrix in the determinant is diagonal in k^, the determinant factors into a 
product of determinants, each with fixed k^. The logarithm of the product is then given 
by an integral over the two-dimensional space of k^ . Since the integrand is invariant under 
rotations in k^, we can write this integral in polar coordinates as 
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where k i 
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After a change of variable p = a/1 + k^/zt^ we obtain the Lifshitz formula for the energy, 
heL^ 



/oo poo 
K'dK j pdplog [{1 - r^r^e-'^^") (l 
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C. Two cylinders 

We now rederive the Casimir energy for two perfectly conducting, infinitely long cylinders. 
The result for one cylinder inside the other has been presented in Ref. [70] and the result 
for both outside each other was presented in Refs. |69l |82] . 

For scattering from outside cylinder a, we have the scattering solutions 

/Ldk' 

n' 

E(/t, x) = Nj^,^^^(K, x) + y 2^ [■'^fc"!n'('*? ^)^a%^n'M,k^nE + ^kf,n'i^^ ^)^a%'^n' E,k^nE] 5 

n' 

(VI.15) 
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FIG. 6: (Color online) Two perfectly conducting infinite cylinders with radii Ra and i?;, are sepa- 
rated by a center-to-center distance d. They can be outside one another, or one may be inside the 
other. 



with boundary conditions E" = and B"*- = on the cylinder surface. L is the length of 
the cylinders, and we are considering the limit L oo. We have 

•'a,k'^n'E,k^nM ~ a,k'^n' M,k^nE ~ ^i 

-r-ee - -"^Slh -A-'U ^ni^aP) 

a,k'n'M,k^nM — L^y'^^ '^zj'^n 



■^Tk'n'E,k,nE — ~'^^{^z ~ ^^^n,\ 



and analogous equations hold for scattering from cylinder h. The energy in Eq. (V.16) is 
given in terms of exterior scattering amplitudes only, 

he r°° 

£ = — y dnhg det (J - J^l'U^'^r^'W^) . (VI.17) 

The matrix inside the determinant is diagonal in k^, so the log-determinant over this index 
turns into an overall integral. A change of variable to polar coordinates converts the integrals 
over K and fc^ to a single integral over p = a/^^ + k^, yielding 

£ = — pdp (log det AT^^ + log detX^) , (VI. 18) 

^ 



where 
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For scattering from inside cylinder a we have the scattering solutions 

E(k, x) = M^"*„(fi;, x) + y 2^ ["'^fcz^n'^'^' ^)-^,k'^n'M,k^nM + ^k'^yi'^^ ^)^,k'^n'E,k 

n' 
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yielding 



''a,k'n'E,k^nM ~ ^ a,k'n' M,k^nE ~ ^ 



(RaP) 



(VI.21) 



^a,k'^n'E,k:,nE — '^L^i^z k^)5n,n' 



In (RaP) ' 

We note that the inside scattering amphtude matrix is the the inverse of the corresponding 
outside result. The energy, expressed in Eq. (V.16), now becomes 
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which contains the appropriate scattering amplitudes for the inside problem, JF** for in- 
terior scattering of object a and JF^'^ for exterior scattering of object b. Using the same 
simplifications as in the outside case, we have 
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where 
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D. Sphere and plate 

In this section we investigate the Casimir interaction of an infinitely thick plate a opposite 
a sphere 6, each with frequency-dependent permittivity and permeability. The geometry is 
depicted in Fig. [7) 

The scattering amplitude for the plate is easy to express in the plane wave basis using 
Eq. (VI. 10). We can apply our result from the plane geometry, Eq. (VI. 12), changing only 
JF^^, which now becomes the scattering amplitude for vector plane wave functions outside a 
sphere. To express the scattering amplitude of the sphere in the spherical vector wave basis, 
we use Eqs. (D.l) and (D.2) and obtain 
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FIG. 7: (Color online) A sphere b of radius R is located opposite a plate a, separated by a center- 
to-surface distance d. 



where P and P' label the polarizations in the planar basis and Q and Q' label the polariza- 
tions in the spherical basis. The normalization factors C\^^p{k) and Cq{k), defined below 
Eq. (B.l) and Eq. (B.7), arise from the definition of the scattering amplitude (see, for exam- 
ple, Eq. (IV.9)). For a sphere with uniform permittivity and permeability, we compute the 



scattering amplitude by solving Eq. (IV.9) in the spherical vector wave basis, which yields 



-r-ee _ nree _ n 
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^blmEI'm'E — ^1,1' ^m,m' 



h{KR)dn{Rii{n,KR)) - fiMRk{KR))ti{n,KR) ' (VI.26) 
ii{KR)dR{Rii{nbKR)) - ebdR{Rii{KR))ii{nbKR) 
ki{KR)dR{Rii{nbK,R)) - ebdR{Rki{KR))ii{nbKR)' 



where Ub is the index of refraction, nbiicn) = \J ebiicn) ^biicK) . The modified spherical Bessel 
functions ii and ki are defined in Appendix |Bj 

= det(I + BA), the energy simplifies to 



Plugging into Eq. (VI. 12) and using det(I 
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where 
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Here k_i_ = |k_|_| and rj, defined in Eq. (VI. 11), is the Fresnel coefficient for reflection of a 



wave with polarization Q. The ratio of Ckj^p(ft) to Cq{k) in Eq. (VI. 25) has opposite signs 
depending on whether P and Q represent the same polarization or the opposite polarization, 
which we have implemented through the term {26p^q — 1). The integration over all angles of 
has already been carried out, which makes A/" diagonal in m and m'. (Although -D;mp,kxQ 
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seems to depend on the angle of k^, the multiplication with its Hermitian conjugate cancels 
this dependence.) 

To leading order for large d/ R, the / = 1 components of the sphere's scattering amplitude 
and the k — limit of the permittivities and permeabilities contribute. The scattering 
amplitude can be expanded to lowest order in terms of the sphere's electric and magnetic 



polarizabilities, T^l^M,imM 
'^R^ and 



and 



ImEAmE 
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\R^ are given in terms of the zero frequency permittivity 



where 
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fib{0)- To leading order the energy is given by 
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can be expressed in terms of elementary functions, but the expressions are too complicated 
to be worth reproducing here. The two functions are plotted in Fig. [8j 

The expression for the energy in Eq. (VI.29) agrees with the results in Ref. [66] for a 
perfect metal plate ea ^ oo and a sphere with general eh- It also agrees with the results in 
Ref. [83] for a perfect metal plate and a perfect metal sphere, oo, and — > cx) . Both 
of these works arrive at similar general expressions for the energy to what we have found 
here; Ref. [HS] combines scattering theory techniques we have used here with the method of 
images, while Ref. [H3] uses Wigner rotation matrices. 

In the calculations of Refs. [66] and [83], when e — > cxd the corresponding /i is set to zero 
to reproduce perfect metal boundary conditions within a low-frequency expansion. (Ref. 



[26] contains the asymptotic Casimir energy formula Eq. (VI.29) in the case where the mag- 
netic permeabilities are set equal to one.) As the plots in Fig. [8] show, however, the perfect 
reflectivity limit of the plate is approached slowly with increasing e^. To compare with ex- 



periments it is thus important to compute the energy (VI. 27) using the actual permittivities 



and permeabilities of the material instead of perfect metal limits. 



E. Cylinder and plate 

As a final example, we investigate the Casimir interaction energy between an infinitely 
thick plate a opposite a cylinder b, each with frequency-dependent permittivity and perme- 
abihty (Fig. [9]). We will focus on presenting asymptotic (d/R — >• oo) results here, but the 
derivation is straightforward to extend to intermediate separations, for which the evaluation 
of the final expression can be performed easily on a computer. We choose the z axis as the 
axis of symmetry of the cylinder and let k_L denote the vector {ky,kz)- 

Just as in the sphere-plate case, it is most convenient to express the scattering amplitude 
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FIG. 8: (Color online) Plots of </>^ (blue, positive) and 0^'^ (red, negative) as functions of 1/eao for 
fixed /Xao = (left) or fixed i^ao = 1 (right). For ^ao = 1 the two functions (j)^ and (j)'^^ approach 
1 rather slowly from the right (perfect metal limit). So, for comparison with experiments, it may 
not be justified to use the perfect metal limit ^ oo of the plate to compute the Casimir energy. 



a 




FIG. 9: (Color online) A cylinder b of radius R is located opposite a plate a, separated by a 
center-to-surface distance d. 



of the cylinder 6 in a plane wave basis by 
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where Ck^p{K) and Cq are defined below Eq. (B.l) and Eq. (B.4), respectively. 

By solving Eq. (IV. 9 ) in a cylindrical wave basis, it is straightforward to find the scattering 
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amplitude of the cylinder, ^t%^nQ,k' n'Q'- For uniform permittivity and permeability, the 
matrix elements are diagonal in and the cylindrical wave index n, but not in TE and TM 
polarization. The expressions are somewhat complicated; since we are presenting asymptotic 
results here, we only need the small-radius expansion, 

H%,nP,k',n' P' = ^ K)^n,n'fk,nPP' + O^R^) , 

fk.OMM = + kl)R^ (1 - fib) , fk^OEE = U^"^ + kDR^ (1 - e;,) , 

Jk,±lMM - 2(l+e6){l+M6) ^ ' Jk.±lEE - 2(1+Ei,)(l+Mb) ^ ' 



f _ -f _ Kk^{tifn,-1) T-,2 r _ f _ r 

Jk^lME — Jk^-IEM — ; JkAEM — Jk^-IME — —JkAME- 
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All other matrix elements (|n| > 1) contribute at higher order in R. It is assumed here that 
ebiicn) is finite. In the infinite conductivity limit (ef, oo) only one of these scattering 
amplitudes contributes; this case is discussed below. 



We next plug into Eq. (VI. 12). As in the case of two cylinders, the matrix inside the 
determinant is diagonal in kz-, so the log-determinant over this index turns into an integral. 
We obtain for the Casimir energy 

£ = — dn dk, logdet (VI.33) 
4vr^ Jo i_oo 



where 



^Jk^nP,n'P' = fk^nPP" / d^y "" ^ ^-= =^ 



(VI.34) 

To find the interaction energy at separations outside of the asymptotic limit, fk^nPP' must 
be replaced by the appropriate scattering amplitude expressions for all n, valid to all orders 



in R. Expanding the logdet in Eq. (VI.33) to first order in A/", we obtain for the interaction 
in the large distance limit d/ R oo, 

^ = + '"^ - ^) - (3 + ^w)x'r'\0, x)] , (VI.35) 

1287ra^ Jo ebo + I 

if the zero- frequency magnetic permeability fibo of the cylinder is set to one. If we do not 
set fibo equal to one, but instead take the perfect reflectivity limit for the plate, we obtain 

^ _ hcLR^ {ebo - /i6o) (9 + ebo + jJ-bo + ebo/^w) fVI 361 

327rc?4 (I + ebo){l + fim) ' ^ ' ' 

Finally, if we let be infinite from the beginning (the perfect metal limit for the cylinder), 
only the n = TM mode of the scattering amplitude, ^b%^oE,k' oe — I7^i^z — K) iog^R/d + 
0{\og~'^{R/d)), contributes at lowest order; the previous expansions of the cylinder's scat- 



tering amplitude in Eq. (VI.32) are not valid. For a plate with zero-frequency permittivity 



Sao and permeability fiao, we obtain for the Casimir energy 
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where 



dx 
l + x 



[rf(0,x)-xrf(0,x)]. 
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In Fig. 10, (j) is plotted as a function of the zero-frequency permittivity of the plate, eao, 



for various zero-frequency permeability values, Hao- 
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FIG. 10: (Color online) Plots of </>^ versus l/cao for fixed values of /iao- The perfect metal limit 
{(p^ = 1) is approached slowly for large jiao^ as in the case of a sphere opposite a plate. For large 



IJLaO the interaction becomes repulsive, which is expected given similar results for two infinite plates 
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APPENDIX A: DERIVATION OF THE MACROSCOPIC FIELD THEORY 



In this appendix we justify the starting point of our derivation, the effective or "macro- 
scopic" field theory in Eq. (II.l ), in order to clarify the causality properties of the permittivity 



32 



e and permeability /i and incorporate the dissipative properties of the materials naturally 
into our formalism. We also show that both non-local and non-isotropic permittivity and 
permeability tensors can be used. 



In place of Eq. (II. 1), we begin with the Lagrangian density of the free electromagnetic 



field plus a coupled system of particles. To be concrete, we imagine that the electromagnetic 
field couples, say, to electrons of charge — e described by the Lagrangian density operator 



E^-B 



2m 



le 



(A.l) 



where ^ and are the fermion field annihilation and creation operators, which are spinor 
functions of space x and time t. We implicitly sum over spins and suppress the spin index 
and we continue to work in = gauge. The electrons' coupling to the lattice, their 
coupling to each other, and their spin-magnetic field interaction can be explicitly taken into 
account by adding, for example. 



Aatticc = -V^"^(t,x)M(x)^/'(t,x) 

-^spin = 7^B(x) ■ ?/;"^(t,x)cr?/;(t,x). 



(A.2) 



2mc 



We imagine that such systems, confined to the various regions of space, represent the objects 
whose Casimir interactions we are calculating. Since the following procedure is quite general, 
we are not limited to systems described by these particular Lagrangians, but modifications 
to our approach may be needed in some situations. 

The electron-lattice and electron-electron interactions are mediated via the quantum 
electrodynamic field, but the relevant wavelengths are substantially shorter than the ones 
dominating the Casimir interaction of different objects, so we can safely approximate the 
short wavelength interactions by effective potentials u{x.) and f (x — x'). 

To compute the total partition function, we exponentiate the time integral of the La- 
grangian and analytically continue the time coordinate t to —ir, yielding 



where 



S^M = \ j dr J rfx (E^ + B^) , 
is the free electromagnetic action and 



■^elcc [A] 



(A.3) 



(A.4) 



(A.5) 



is the electronic partition function. The Euclidean electronic action is 



dr I d^ [ hipdrip H 

J \ 2m 



+ interactionsj , 
(A.6) 
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where the interaction terms are taken from Eq. (A. 2), including at least Cutucc to keep the 



electrons confined to the various objects. Here the creation and annihilation operators ip"^ 
and ip go over to Grassman path integral variables ip and ip. 

Next, we expand the partition function of the electronic system, 



where the second derivative of Zgiec, 
Kij{T - r',x,x') 



,,,r-r',x,x')A^(r,x)A^(r',x')j , 

(A.7) 
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Zeiec[0] M*(r,x)Mi(r',x') 
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A=0 



only depends on the difference in imaginary time r — r'. The linear term has been omitted 
in this expansion because it vanishes for systems with no mean currents. We then obtain 
the well-known Kubo formula for electrical conductivity, [85] 



Cr,^(zC/t„,X,x') 
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The T superscript indicates that this is the imaginary-time ordered response function. The 
retarded response function can be obtained by the substitution icKn — > a; + iO"*". The 
conductivity is related to the permittivity by 



X, X') = CKn (e^(iCK„, X, X') - (5ij(5(x - x')) . 



(A.IO) 



After substituting into Eq. (A.7), we obtain 



Zelec ~ ^clcc[0] I 1 + 2 X] / C?X(ix'E* (ic/tn, x) ■ (M(x - x') - {iCKn, X, x')) E(iCKn, x') 
\ n=— 00 



(A.ll) 



and, finally, after reexponentiating and plugging into Eq. (A. 3) we obtain 



~ Zciec[0] / PAexp —P ( / (ixcix' E*(icK„, x) ■ e"^(ic/€„, X, x')E(2CK„, x') 



(ix B* (icKr, , x) ■ B (icK,„ , x"! 



(A.12) 



The imaginary-time ordered response function can be obtained from the retarded real time 
response function, which is experimentally accessible, by analytic continuation to imaginary 
frequencies. The relationship e^(icfi;„, x, x') = e^(zc|K„|, x, x') between the two response 
functions and the symmetry of the retarded response function in the indices (z, x) and (j, x') 
is discussed in Refs. |86l |87] . 

To be able to consider two objects as separate and distinct, we assume that e^{icKn, x, x') 
is zero when x and x' are on different objects, separated by the vacuum or a medium. This 
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is justified even for small gaps between the objects, since tunneling probabilities decay 
exponentially in the gap distance. 

Although the dissipative part of the electric response (the conductivity) appears in the 



abov e equations, there is not actual dissipation at zero temperature. Rather, Eqs. (A.7)- 
(A.9) merely show that the conductivity cr'^, which can be measured in an experiment. 



can be related to the fluctuations that exist in a system in the absence of an external 
perturbation. The size of these fluctuations is represented by the second derivative of the 



partition function, Eq. (A. 8) 



Although the magnetic susceptibility is usually negligible compared to the electric suscep- 
tibility, we have kept the permeability function fi in our derivations of Casimir interactions. 
If the spin-magnetic field coupling term £spin replaces or dominates the standard coupling 



between charge and electromagnetic field in Eq. (A.l), an analogous procedure introduces 



the magnetic permeability function into the action. In that case, the partition function of 



the matter field, the analogue of Eq. (A.7), has to be expanded in B* = (V x A)* instead 
of A\ 



APPENDIX B: GREEN'S FUNCTION EXPANSIONS AND MODIFIED 

EIGENFUNCTIONS 

In this section we supply the Green's function expansions for imaginary frequency icK 
and the associated modified eigenfunctions in various bases. After analytically continuing 



the wave functions in Eq. (III. 8), it is convenient to redefine them by multiplication by an 
overall factor in order to obtain the conventional definitions of the modified wave functions. 
Since E''^*(k, x) = E°"*(fi:, x), it suffices to supply the modified regular and outgoing wave 
functions. 

Electromagnetic vector waves are typically divided into TE and TM modes. It is custom- 
ary to name the TE wave functions M(ci;,x) and the TM waves N(ci;,x). Often TE modes 
are referred to as magnetic modes, labeled by M, and TM modes are referred to as electric 
modes, labeled by E. 



1. Green's function and eigenfunctions — plane wave basis 

We choose the z axis as a convenient symmetry axis and let k_L denote momentum 
perpendicular to this axis. The free dyadic Green's function is 

Go{icK,, X, x') = 

dk^ f Ck,M(«:)M°'f(«:, x) ® Mj^f (/s:, x') + C^^e{k)NI-^{k, x) ® Nj^f («:, x') if ^ > ^ 

(B.l) 
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Here, Ckj_M{K) — — , \ ^ — (—l)Ckj^ and the modified vector plane wave functions 
are given by 

ML'!(/.,x) = X C!(/^,x)z, M-;(/.,x) = X C(/.,x)z, 

Kli^, x) = ^ V X V X («:, x)z, N-*(«:, x) = ^ V x V x x)z, 

(B.2) 

in terms of the modified scalar plane wave functions, 



As discussed in the text, the labels "reg," "out," and "in" are not really appropriate for 
plane wave functions, but the mathematical results all carry over. 



2. Green's function and eigenfunctions — cylindrical wave beisis 



Again, let us take the z axis as a convenient symmetry axis and let p = -\--xp' be the 
distance to the z axis. The free dyadic Green's function is given by 



[ro(«CK, X, X 

CmM-:,(«:, x) ® M-^;(«:, x') + C^N-* («:, x) ® N-*;*(k, x') if p > p' 
CmM^^J/., x) Mjr;(«:, x') + Ce^':%{k, x) ® N'^;„(«:, x') if p < p' 

(B.4) 



n ^ 



(B.5) 



Here Ce — ^ — (— 1)Cm, the vector cylindrical wave functions are given by 

M?>.x) = ^;^^Vx0->,x)z, 

MS(k,x) = J V X r^nMi, 

N°,:U«^,x) = X V X Cn(«^,x)z, 

and the modified cylindrical wave functions are 

Cn(«^, x) = 7„ (paAFT^) e^'=^^+^"^ Cn(«^, x) = X„ (pv/^FT^) e^'=^^+^"^ (B.6) 

where J„ is the modified Bessel function of the first kind and Kn is the modified Bessel 
function of the third kind. 
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3. Green's function and eigenfunctions — spherical wave basis 

In spherical coordinates the free dyadic Green's function is given by 



uo{ichi, X, x') 



E 

Im 



CMin)MZ'M ® M[r(«:,x') + Ce{^)NZ'M ® KSi^^^') if |x| > |x'| 
CM(«:)M;:f(«:,x) ® M|-(«:,x') + C^(«)N[;^(«:,x) ®Ni-(«:,x') if |x| < |x'| 



(B.7) 

Here Cm{,k) = k = {—1)Ce{i^), the vector spherical wave functions are 

= Tim^ C'('^,x)x, MZ'M = ^;^V X 0r('^,x)x, 
^^-^^'^^ = ^ C'(«^,x)x, NZ'M = ^^^^ X V X 0r(«^,x)x, 

(B.8) 

and the modified spherical wave functions are 

cPTr^M = ii{k\^\)Y,U^), ^ZnM = A:^'^ |x| (x) , (B.9) 



where ii{z) = a/^-^/+i/2(-2) is the modified spherical Bessel function of the first kind, and 
ki{z) = \/ ^Ki^i/2{z) is the modified spherical Bessel function of the third kind. 



4. Green's function — elliptic cylindrical basis 



In order to study geometry and orientation dependence of Casimir interactions, it is 
helpful to be able to study objects with reduced symmetry. In Ref. [88], this formalism was 
applied to spheroids in scalar field theory. Unfortunately, the vector Helmholtz equation is 
not separable in spheroidal coordinates as it is in spherical coordinates. While the analogous 
vector spheroidal harmonics can still be constructed, the scattering matrix for a perfectly 
conducting spheroid is not diagonal, although it can be obtained from a more elaborate 
calculation |89]. For a perfectly conducting elliptic cylinder, however, the vector scattering 
problem is separable, so we describe that case here. Throughout this section, we use the same 
normalization and conventions as in Ref. [90], in which all functions in elliptic cylindrical 
coordinates have the same normalization as their circular analogs. As a result, all the 
functions inherit the usual completeness and orthonormality relations and approach their 
circular analogs in the limit of long wavelength. 

In elliptic cylindrical coordinates, the z coordinate is unchanged, while the components 
of xx become x = a cosh.fi cos 6 and y = asinh/xsin6', where the interfocal separation of 
the ellipse is 2a. Far away, 6 approaches the ordinary angle in cylindrical coordinates and 
|x_l| ^ le'^. Separation of variables in these coordinates yields angular and radial Mathieu 
functions for 9 and /i, respectively. The even and odd angular Mathieu functions are cen{0, 7) 
with n > and sen{d, 7) with n > 0, which are the analogs of cos n6 and sin n6 in the circular 
case. (We used a complex exponential basis for the circular case, but it could equally well be 
represented in terms of sines and cosines.) The angular functions now depend on the wave 
number through the combination 7 = —{k"^ + n'^)a? /2. The corresponding radial functions 
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are now different for the even and odd cases and depend on 7 and separately rather than 
through a single product of the two. The even and odd modified radial Mathieu functions 
of the first kind are denoted Iem{fJ','j) and l0m{l^,l) respectively, and the even and odd 
modified radial Mathieu functions of the third kind are denoted KemifJ^yj) Komi^ijj), 
respectively. 

We then obtain the same results as in cylindrical coordinates, but now with 

Cne(«:. x) = /e„(/i, 7)ce„(^, 7)e''^^^ €Tnei^, x) = Kenifi, l)cen{9, 7)e''=^^ 
Cno('^,x) = /o„(/x,7)3e„(^,7)e^'^^ <Pi:U^,^) = Ko4f^,^)se^{e,^)e''^\ 

For numerical calculation the required Mathieu functions can be efficiently computed using 
the C++ package of Alhargan [91] 192] . Analogous replacements convert the translation 
matrices and wave conversion matrices described below into this basis. 



APPENDIX C: TRANSLATION MATRICES 



In the following, we list the translation matrices that make up X*-^, defined in Eq. (III. 14). 
The definition of the vector Xjj, which points from the origin of object i to the origin of 
object j, is illustrated in Fig. [TJ 



1. Plane wave basis 

Plane waves are eigenfunctions of the translation operator, which does not mix TE and 
TM vector plane wave functions. 

If the z coordinates of object i are smaller than those of object j, then —V*-' is the only 
nonzero entry in X*-'. Taking Xjj to point from the origin of object i, Oi, to the origin of 
object j, Oj (that is, upward), we obtain 



Cp,kYP' = e-^''-^-'^-v/^^-^5(2)(k^ - kl)5p,p,. (C.l) 



If i is located above j, then — W-'* is the only nonzero entry in X*-'. The vector Xjj points 
upward from Oj to Oi, and we have 

■\Aj3i - VJ" C'k|P(«:) 

' ^^1, , p pi ~ "^t' pi 



L2 

Since the matrix is diagonal in k_|_ and P, the factor ^!^^^^f^-. cancels. 



2. Cylindrical wave basis 

Translations do not mix the TE and TM modes of vector cylindrical wave functions. 
They are constructed by taking the scalar cylindrical wave function, multiplying by the 
unit vector z, and performing one or two curl operations. A TE vector cylindrical wave 
function is perpendicular to z, while the curl of a TM vector cylindrical wave function is 
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perpendicular to z. Expanding any of the two vector wave functions around any other point 
in space must preserve its orthogonahty property with respect to the constant vector z. So, 
the two arc not mixed by the translation matrix. 

If two objects i and j are outside of one another, —W'' is the only nonzero submatrix of 
X'-'. Again, let X^j point from Oj to Cj. We have 

KnP,Kn'P' = Kn-n' {p^r^A^W^) e-^^^'^^-^^"-"')^^'' (- l)"'5p,p. f 5(A;. - K), (C.3) 

where |Xjj _|_| is the distance of X^j to the z axis, i.e. the length of the projection onto the 
x-y plane, and 6ji is the angle of X^j in the x-y plane. 

When object i is enclosed inside the surface of an infinite cylinder, inside object j, sub- 
matrix — V'-^ is the only nonzero entry in X*-' . We have 

where Xj^ points from Oi to Oj, 

If the roles of i and j are reversed, then —W^'' is the nonzero submatrix of X*^, with 

'^^k^nP,k'n'P' ~ ^ k'n' P' ,k;,nP ri 

' ' <^P' (C.5) 

Since the matrix is diagonal in P, the factor ^ cancels. 



3. Spherical wave basis 



The TE vector wave functions are orthogonal to the radius vector x. Since the same 
vector wave function cannot also be orthogonal everywhere to the radius vector of a shifted 
coordinate system, TE and TM polarizations mix under translation. 

Suppose object i and its origin are outside a spherical separating surface, which encloses 
j. The nonzero submatrix of X*-' is —W\ with 



i^Lmmm - i-^r^' E [^(^ + 1) + ni' + 1) - I'V" + 1)] J "^'^lafiw^^^ 



X 



I" 

I /' l"\ f I V I" 
0/ \m —m' m! — m 



l{l + !)/'(/' + 1) 

fc/// ( K I Xj j I ) yj//^-^/ (Xji ) , 



^l'm'E,l: 



+ !)/'(/' + 1) 



^2 {Km^Vm'lm+lO'^ji) + \m'^l'm'lm-lP^ji)) 



TP'' 

^l'm'M,lmE 



^l'm'E,lmM^ 



'E,lmE ^H'm'M,lmM^ 



where 



Avm'lmi^^ji) - (-1)-+' J2 VM'^l + 1)(2/' + 1)(2/" + 1 

I" 

'I r I" 



X 



/ r 



I" 



0/ \m —m m—m 



kill (kIXji I )l/"TO-m' (X 



■31') 



(C.6) 
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and = ^{lTm){l±m + l). 

The translations between regular waves are described by the matrix elements 

viLmmm - (-1)- E m + 1) + ni' + 1) - ni" + 1)] 

" (o o) {L -m' J^H«|X.|)(-l)^'V.„KX.), 



^2 {^im^l'm'lm+lO^ij) + \rn^l'm'lm-l0^ij)) 



+ {^im^l'm'lm+lO^ij) — \m^l'm'lm-l0^ij)) + ^ Bif^'irnO^ij] 

^l'm'M,lmE ~ ^ V m' E ,lmM ^ ^l'm'E,lmE ~ ^ V m' M ,lmM ^ 

where 

B.m'imi^a) = i-^r E V47r(2/ + 1)(2/' + l){2l" + 1) 



^ fn n n) ^ ) v('^|X.,|)(-l)'"ll"„'(X,,) (C.7) 

\0 / \m — m m —mj 



and = =f m)(/ ± m + 1). 
The matrix W-^* is related to V-^*, 

V^*, of course, is the same as V*-' with Xj^ replaced by X^j. To be more specific, the 
elements correspond in the following way, 

'^^l'm'M,lmM ~ ^lmM,l'm'M ^ ^ ^ V m' E ,lmM ~ ImM ,1' m' E ^ (C Q\ 

ywJi _ _yji* lyji _ yJi* ^ ' ' 

^^Vm'M,lmE— ^lmE,Vm'Mi ^ ^ V m' E ,1'mE — ^ ImE ,V m' E ■ 

APPENDIX D: WAVE CONVERSION MATRICES 

It is not necessary to express all the objects' scattering amplitudes in the same basis. 
Here, we supply the matrices that convert modified vector plane wave functions to spherical 
or cylindrical vector wave functions. 

1. Vector plane wave functions to spherical vector wave functions 



'™ (D.l) 

Im 
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The conversion matrices are obtained from the decomposition of a plane wave in spherical 
coordinates, 



47r(2/ + l)(/-m)! |k^| - / ATT^/ \ 

DlmE,-k.j_E = DimM,kj_M, DimM,kj_E = —DimE,k±M, 

where P™ is the associated Legendre polynomial and prime indicates the derivative of PI" 
with respect to its argument. 



2. Vector plane wave functions to cylindrical vector wave functions 

The cylindrical vector wave functions are defined as before, but now we consider regular 
vector plane wave functions that decay along the — x axis instead of the — z axis, 

^ \ (D.3) 

The vector plane wave functions can be decomposed in vector cylindrical wave functions, 

n 
n 

using the conversion matrix elements 



(D.4) 



Dk^nE,k±E — Dk^nM,k±M , D};^nM,k±E — —Dk^nE,k±M, 

where ^ = ^^^^-^ and k_L = {ky, k^). 
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